Anisotropic sub-Doppler laser cooling in dysprosium magneto-optical traps 
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Magneto-optical traps (MOTs) of Er and Dy have recently been shown to exhibit populationwide 
sub-Doppler cooling due to their near degeneracy of excited- and ground-state Lande g factors. 
We discuss here an additional, unusual intra-MOT sub-Doppler cooling mechanism that appears 
when the total Dy MOT cooling laser intensity and magnetic quadrupole gradient increase beyond 
critical values. Specifically, anisotropically sub-Doppler-cooled cores appear, and their orientation 
with respect to the quadrupole axis flips at a critical ratio of the MOT laser intensity along the 
quadrupole axis versus that in the plane of symmetry. This phenomenon can be traced to a loss 
of the velocity-selective resonance at zero velocity in the cooling force along directions in which 
the atomic polarization is oriented by the quadrupole fleld. We present data characterizing this 
anisotropic laser cooling phenomenon and discuss a qualitative model for its origin based on the 
extraordinarily large Dy magnetic moment and Dy's near degenerate g factors. 

PACS numbers: 37.10.De, 37.10.Gh, 37.10.Vz 



I. INTRODUCTION 

Ultracold dipolar gases are emerging as an important 
tool for exploring strongly correlated physics in a system 
with easily and rapidly tunable parameters. Highly mag- 
netic atoms provide the ability to realize strongly corre- 
lated systems in the absence of polarizing fields, which 
will be important for exploring, e.g., quantum liquid 
crystal physics [1 . Additionally, degenerate fermionic 
and bosonic dipolar gases are expected to exhibit a 
wide range of interesting quantum phases and dynam- 
ical properties [2]. Dysprosium, the most magnetic laser 
cooled and trapped gas [3J |1| , is especially exciting since 
its /i = 10 ^B magnetic dipole moment — coupled with 
its large mass and abundance of bosonic and fermionic 
isotopes — is sufficiently strong to allow wide exploration 
of ultracold dipolar physics [3H5] . 

Sub-Doppler cooling [S] in an optical molasses [7] sub- 
sequent to MOT capture is a crucial step toward creating 
degenerate gases in, e.g., alkali atomic systems such as 
Rb and Cs. Sub-Doppler mechanisms fail to cool the ma- 
jority of atoms [F inside these MOTs because the ground 
state Qg factor is substantially different from the excited 
state's (for ^'^Rb, Ag/g^ = 34%) 0; population- wide cr+- 
a~ sub-Doppler cooling breaks down in the longitudinal 
magnetic fields of a typical MOT. Notable exceptions [10] 
are the highly magnetic lanthanides Er (7 /i^) and Dy, 
whose repumper-less MOTs [Sj [TT] operate on optical 
transitions whose ground and excited state g factors are 
the same to within 0.3% (1.7%) for Er (Dy) [12]. At 
low MOT beam intensity, both the Er and Dy MOTs ex- 
hibit in-situ, population-wide sub-Doppler cooling |Hll3j. 
The minimum temperatures >200 fiK are, however, more 
than 10 times larger than typical Rb and Cs optical mo- 
lasses sub-Doppler cooled temperatures, but are similar 
to sub-Doppler temperatures in the ®^Sr MOT [Hj . 

We describe here a novel sub-Doppler cooling mech- 
anism in the highly magnetic Dy MOT system that we 
first reported in Refs. [2111]. We expand on these observa- 
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FIG. 1. (color online). '^^''Dy energy level structure 12 . The 
MOT and Zeeman slower employ the strongest laser cooling 
transition between the even parity (red) ground state and the 
odd (black) excited state at 421 nm. For this transition, the 
g factor of the ground (excited) state is 1.24 (1.22) and the 
Doppler cooling limit is T£ = 768 fiK. J is the total electronic 
angular momentum. 



tions with additional measurements and offer a plausible 
explanation for this phenomenon based on the theory of 
velocity selective resonances (VSR) in a highly magnetic 
gas [T51 - fT7] . While this mechanism only cools a small 
fraction of the atoms to ^10 /iK ultracold temperatures, 
its existence is a first example of novel behavior arising 
from laser cooling the most magnetic atom. 



II. DY MOT REGIMES 

We choose to describe here the Dy MOT system in 
which all pairs of MOT beams are well-aligned in mutual 
orthogonality and retroflection [TS'. All data are taken 




FIG. 2. (color online). (a,c,e) Time-of-flight expansions and 
(b,d,f) in-situ images of the three striped MOT regimes ex- 
hibiting anisotropically sub-Doppler cooled cores. (a,b) Sym- 
metric/cross MOT with Iz/Ip = 1.4. Notice the faint vertical 
and horizontal stripes crossing at the MOT center. (c,d) Pro- 
late striped MOT with IJIp = 0.5. (e,f) Oblate striped MOT 
with h/Ip = 2. 



with the experimental apparatus for Dy MOT produc- 
tion described in Ref. ||j. Briefly, a high-temperature 
oven along with transverse laser cooling produce a colli- 
mated atomic beam. A Zeeman slower, operating with 
1 W of laser power on the 421-nni cycling transition, de- 
celerates the beam before MOT capture on the 421-nm 
line (see Fig. nl). No repumper is necessary since the 
highly magnetic atoms remain confined in the MOT's 
magnetic quadrupole trap as population decays through 
metastable states [3l|4j. We focus here on the bosonic iso- 
tope ^® Dy with zero nuclear spin /, but similar results 
were observed in the fermionic ^^-^Dy MOT (/ = 5/2). 

The Dy MOT can be formed in two classes of oper- 
ation, striped and stripeless, and the striped MOT can 
be further classified into three regimes. The stripeless 
Dy MOT is similar to the Er MOT reported in Ref. p] . 
in that the majority of atoms are sub-Doppler cooled 
to ^200 /iK. These MOTs are obtained at low inten- 
sity / = It/Is ^ 0.17, where It = Iz + 2Ip is the total 
MOT intensity (one beam pair in z and two in p) and 
Is is the saturation intensity in the MOT; Ig — 2.7 x 58 
mW/cm^ g]. Moreover, the stripeless ^'^^Dy MOT forms 
at low MOT magnetic quadrupole field VzB < 12 G/cm, 
where V zB k, 2V pB and z and p are the directions along 
the quadrupole's axis and plane of symmetry, respec- 
tively. 

The temperature and density distributions in the three 
striped MOT regimes differ greatly from the stripeless 
MOT in time-of-flight. The temperature distribution is 
anisotropic at large / and V^i?, and Fig.[2]shows charac- 
teristic time-of-flight expansions of the striped Dy MOT. 
In these images, MOT beams and the quadrupole field 
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FIG. 3. (color online). Temperature characterization of the 
prolate striped Dy MOT [Fig. [2fc)] versus normalized in- 
tensity I and detuning 5. The orange squares (blue dia- 
monds) are the Tmajority along p [z). Light blue curves show 
the Doppler cooling limit [7j for the MOTs' parameters: (a) 
5 = -1.2, V^B = 20 G/cm; (b) / = 0.2, V;,B = 20 G/cm. 



are extinguished at td = 0, and the in-situ images are 
taken at t^ < 0. We observe in time-of-flight a low- 
population core surrounded by a more populous and hot- 
ter outer shell, and we designate the group of outer shell 
atoms as the majority while the inner core as the mi- 
nority. Parametrically driving the MOT could reveal the 
minority component's in-situ size J19j . 

Typical high-population MOTs are formed with / = 
0.2 and i5 = A/T = —1.2, where A is the detuning from 
the r = 27r • 32 MHz, 421-nm transition. These typical 
MOT parameters correspond to Td's of approximately 
1 mK. At large MOT intensity, the temperature of the 
majority component ^majority is consistent with T^i, as 
shown in the p (orange square) data of Fig. [s] and as 
discussed in detail in Ref. [1] . The Tmajority data in p are 
derived from the broader of the two Gaussians employed 
to fit intensity versus p [see Fig. [lib)]. Temperature data 
for z (blue diamond) are derived from single Gaussian fits 
to intensity versus z [see Fig. |4Fc)]. 

The larger z temperature data are a result of the convo- 
lution of the Doppler-cooled majority temperatures with 
the min ority velocity distribution along z (described in 
Sec.|Vl|). We note that Tmajority data in z below / « 0.17 



are roughly equal to the p data and are not shown. This 
low-intensity, sub-Doppler cooled region is also the re- 
gion in which stripes cease to appear in time-of-flight ex- 
pansions: In the Dy system, the stripeless regime is the 
low-intensity limit of the striped regime. We note that 
the sub-Doppler-cooled Er MOT, which did not exhibit 
stripes, was also operated below I — 0.17 in the data of 

Ref. nni. 



III. QUANTITATIVE CHARACTERIZATION 
OF ANISOTROPIC SUB-DOPPLER COOLING 

The minority component assumes either a prolate, 
symmetric/cross, or oblate shape in time-of-flight de- 
pending on whether the ratio of intensity in the MOT 
beams Iz/Ip is <1.4, ~1.4, or >1.4 (see Fig. [2]). MOT 
populations are largest in the prolate, Iz/Ip < 1-4 



td = 7 ms 



lz/lp = 0.5 




FIG. 4. (color online). Prolate striped MOT time-of-flight 
expansion at i^ = 7 ms. (a) Double-Gaussian fit, orange line, 
along a z intensity integration, (b) Intensity integration along 
p is consistent with a single Gaussian fit. 



regime, whose minority component temperature Tminority 
and ratio of minority atoms to total MOT population 

A^minority/iVtotal IS showu in Fig. [S] [20] . 

^minority and A^minority are extracted from a double- 
Gaussian fit to time-of-fight data such as that presented 
in Fig. [4] The temperature of the minority component 
is anisoptropic, leading to the vertical stripe structure in 
Figs.l2jc) andllFa). The temperature in z is hotter than 
the majority atoms [see the z (blue diamond) data of 
Fig. p] , while the temperature in p is sub-Doppler cooled 
to temperatures as low as 10 /iK, a factor of 10 less than 
the sub-Doppler temperatures observed in the stripeless 
regime. The p temperature of the minority component in 
the prolate MOT is well described by a ID characteristic 
sub-Doppler scaling law O [TJ] : 



T„=Tn 



a 



hvi 

2kB\5\ 



(1) 



where [C„+„-] = [0.1(0.09), 0.55(0.20)] in Figs, [ga) and 
(b), respectively. These values of [C^+cr-] are consis- 
tent with those of the Er MOT [T3], but smaller than 
in the ^"^Sr MOT [M]. However, Tq is <10 pK in both 
fits and is ~10x lower than in the Er, Sr, and stripeless 
Dy MOTs. This 2D, anisotropic sub-Doppler cooling — in 
the case of the oblate striped MOT, ID j21] — is a much 
more effective cooling mechanism than the one found in 
the stripeless sub-Doppler-cooled Er and Dy MOTs. 



IV. TRANSITION BETWEEN PROLATE AND 
OBLATE STRIPED MOTS 

We explore the transition between the prolate and 
oblate striped Dy MOT in Fig. [6] As the MOT's Ijlp 
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FIG. 5. (color online). Tminority in p of the prolate striped 
MOT h/If, = 0.5 [Fig. [2];c)]. Blue curves in (a, b) are fits 
of the data to the sub-Doppler cooling expression Eq. [l] (d- 
f) Fractional population of the minority component, (a, d) 
S = -1.2, V^B = 20 G/cm; (b, e) 7 = 0.2, V^B = 20 G/cm; 
(c, f) 1 = 0.2, S = -l. 



ratio is tuned from 0.5 to 2, we note the occurrence — 
defined by visibility of the double-Gaussian structure 
above image noise — of stripes in a series of 10 realiza- 
tions. The fraction of these images exhibiting a stripe 
structure is noted in Fig. |6] Below Iz/Ip — 1, the pro- 
late stripe is always observed, while the oblate stripe is 
always observed above Iz/Ip = 2. There is a smooth 
transition between the regimes in which successive MOT 
realizations may exhibit either a prolate or oblate stripe, 
and we can offer no explanation for either the smoothness 
of the transitions or their slight asymmetry. 

Between these two regimes there is a third: We most 
often observe a spherically symmetric minority compo- 
nent in the range of Iz/Ip from 1.3 to 1.5. This is the 
regime of the power balanced and aligned MOT discussed 
in Ref. j3] and is less populous, though more dense, than 
the prolate regime. A prolate-oblate cross is observed 
[see Fig. l2ja)] in the region around the critical ratio 
Ir = Iz/Ip = 1-4, though not reproducibly, and the 
spherically symmetric core is more prevalent. 

No stripes are observed below the critical intensity 
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FIG. 6. (color online). Visibility of stripe regimes — prolate, 
symmetric/cross, and oblate — versus the ratio of MOT beam 
power Iz/Ip- Data are taken for I — 0.22(0.01), S = —1, and 
VzB = 20 G/cm. All stripes vanish below Ic = 0.17. Each 
point is the average of 10 MOT realizations with td = 6.5 ms. 



Ic < 0.17, which coincides with the appearance of ma- 
jority component (population-wide) sub-Doppler cooling 
[see Fig.ls] . The prolate stripe is observed to vanish below 
a critical VzBc ~ 12 G/cm, as shown in Fig. IsFc). The 
minority component temperature and stripe population 
fraction rise below V^Bc — which may be interpreted as 
a blending of the minority and majority components — 
until the MOT is firmly in the stripeless, isotropically 
sub-Doppler-cooled regime. 



V. VELOCITY SELECTIVE RESONANCE 

PICTURE OF ANISOTROPIC SUB-DOPPLER 

COOLING 

This section presents a qualitative explanation of the 
anisotropic sub-Doppler cooling mechanism. Two prop- 
erties of Dy are crucial to this explication: its large 
magnetic moment and the near equal ground and ex- 
cited state g factors {Ag « 0) of the cooling transi- 
tion. The latter is important for the following reason. 
In standard ct^-ct^ sub-Doppler cooling, the linear op- 
tical polarization serves as the quantization axis since 
the optical pumping rate jp is typically greater than 
the Larmor precession rate ujl = IJ-V Br/h, where r is 
the distance of the atom from the magnetic quadrupole's 
center. In low magnetic fields, ground state population 
imbalance due to non-adiabatic following of the quanti- 
zation axis induces the differential scattering of ct- light, 
which then leads to sub-Doppler cooling to zero mean 
velocity [51 |Tj. However, in large longitudinal mag- 
netic fields, a non-zero Ag results in an "unlocking" of 
the Doppler and sub-Doppler-cooling mechanisms, and 
sample temperatures rise to the Doppler-cooling limit. 
(Hence the failure of population-wide sub-Doppler cool- 
ing in, e.g., Rb MOTs and the need for field nulling to 
the mG level when sub-Doppler cooling in Rb optical 
molasseses.) By contrast, Er and Dy MOTs remain sub- 
Doppler-cooled [H [T3j despite the large intra-MOT mag- 



netic fields because Ag w in both. This explains the 
existence of intra-MOT sub-Doppler cooling in Er and 
Dy, but not the origin of the anisotropic regime. 

The anisotropic sub-Doppler-cooling mechanism may 
be qualitatively understood in the velocity-selective reso- 
nance (VSR) picture of ID sub-Doppler cooling [T5] when 
augmented to account for large magnetic fields, or equiv- 
alently, large magnetic moments |22| . In this picture, 
most types of sub-Doppler cooling may be understood as 
arising from the momentum transfer by a coherent two- 
photon process between ground state sub-levels (labeled 
by m) 15J. These coherent Raman transitions can oc- 
cur when the difference in Doppler-shifted frequencies of 
counterpropagating light beams seen by an atom equals 
the Zeeman shift between levels separated by Am — n. 
Specifically, VSRs occur when Sk ■ v — njisB, where v 
is the atom's velocity and 5k = 2fc is the difference in 
the wavevectors of the cr^-tr" light. The allowed n's de- 
pend on the relative orientation of optical and magnetic 
field vectors, and it is this selection rule that is the origin 
of the anisotropic sub-Doppler cooling in the large field 
(large magnetic moment) regime. 

To see how this selection rule leads to anisotropic cool- 
ing, we first note that the atoms' polarization aligns along 
the local magnetic field in the presence of a large mag- 
netic field (or large moment), defined as 7p <C wl ^ F; in 
this large Zeeman energy regime, the magnetic field now 
serves as the quantization axis, not the optical polariza- 
tion |16j . Now we examine the case in which the mag- 
netic field is aligned parallel to the fc's of a pair of cr+-cr~ 
beams. Along the spatial directions defined by the fc's of 
the cooling light, the zero-mean- velocity VSR in the cool- 
ing force disappears {n ^ 0) even though finite-velocity 
VSRs remain {n>2 are allowed). This is because there 
is no component of the polarization along the magnetic 
field (azimuthal symmetry about k is preserved), and 
thus no TT transitions are allowed, which are necessary 
for n = and n = 1 VSRs. In this case, sub-Doppler 
cooling to zero velocity along directions close to k is not 
possible when B \\ k. 

However, not all sub-Doppler cooling is suppressed. A 
zero-velocity VSR remains in the cooling force along di- 
rections in which the local magnetic field is perpendicular 
to the cooling light's fc's (for B _L fc all n are allowed). 
This spatial modulation of allowable VSRs (n's) is the 
origin of the anisotropic nature of sub-Doppler cooling 
observed in the highly magnetic Dy system. 

In summary, a^-a^ sub-Doppler cooling is preserved 
(fails) along directions perpendicular (parallel) to the lo- 
cal magnetic field in systems with large magnetic fields 
or large magnetic moments. Figure mjb) illustrates this 
idea for the prolate striped MOT: sub-Doppler cooling 
occurs along p due to the green fc^ a'^-a~ beams, but 
not along z due to atomic polarization orientation par- 
allel to the magnetic field. The red fe^ cr^-cr" beams 
only induce VSRs with non-zero velocity (e.g., n > 2) in 
the cooling force, and these beams cause the atoms to 
accelerate along z due to the magnetic quadrupole field 




FIG. 7. (color online), (a) Orientation of the MOT's magnetic 
quadrupole field, (b) Illustration of the anisotropic cooling 
in the prolate MOT including the MOT magnetic (B) field 
orientation, cooling beam wavevectors fc, and polarizations 
a near the z-axis. 



gradient. The analogous mechanism holds for the oblate 
MOT when the magnetic field (quantization axis) is along 
p instead. 

Prolate and oblate stripes oriented with respect to 
the quadrupole axis dominate because the only regions 
wherein the magnetic field is parallel to the fc-vectors are 
those in the plane or along the axis of symmetry of the 
quadrupole field. The regions' small relative volume is 
the origin of iVminorityZ-^totai < 10%. We note that this 
a~^-a~ mechanism is distinct from a^-a^ magnetically 
induced laser cooling [7]. 



VI. ANISOTROPIC SUB-DOPPLER COOLING 

IN THE DY SYSTEM 

Although this description does not predict the prolate- 
to-oblate transition, our quantitative measurements in 
the prolate regime are consistent with this picture. 
Specifically, when Iz/Ip < I^, I > L, and V^B > V^Bc 
such that a prolate striped MOT forms, we observe 
higher-than-Doppler-limited temperatures along the z di- 
rection, which lies parallel to the local magnetic field 
(Fig. w|. As mentioned in the previous section, this is 
because along the B \\ k direction the atoms accelerate 
as they move away from the trap center due to ever in- 
creasing magnetic fields and VSR velocities. An estimate 
of these resonant velocities — based on ii = uj^/k [15 in 
the 20 G/cm gradient field along z — are > Tjj, and are 
thus consistent with the absence of any sub-Doppler cool- 



ing mechanism along this direction. Specifically, outside 
a core of radius 30 fim — roughly the zone outside which 
anisotropic sub-Doppler cooling occurs; see below — the 
velocities are >0.35 m/s, resulting in ID temperatures of 
>2.5 mK. Such temperatures are larger than the To « 1 
mK limit for the typical Dy MOT parameters (see blue 
curves in Fig. pi) . 

Anisotropic sub-Doppler cooling is possible for the Dy 
MOT on the 421-nm, Ag « transition because 7p <C 
wl <C r is satisfied throughout the majority of the Dy 
MOT's quadrupole field. Specifically, Dy's optical pump- 
ing rate jp ~ /r/(l + 1 + 46^) is approximately equal to 
the light shift frequency, and the inequality is satisfied 
with the range of values [7^, wl, F] — 2tt ■ [0.8, 1 — 28, 32] 
MHz present in a typical Dy MOT with VB = 20 G/cm 
and 0.04 < r < 1 mm. In addition, "fp/uj^ < 1 for 
VB > 10 G/cm, which is close to the measured critical 
value Wz^c/'2 — VpBc ~ 12 G/cm. We do not have 
a simple estimate for Ic ~ 0.17 and I^ = 1.4, but we 
expect that 2D numerical simulations of the Dy MOT 
could explicate the origin of these specific values as well 
as the underlying mechanism behind the stripe orienta- 
tion transitions. 



VII. CONCLUSION 

In conclusion, we note the possibility that other 
highly magnetic lanthanides could exhibit anisotropic 
sub-Doppler cooling for large enough cooling beam in- 
tensity and field gradient: This novel cooling mecha- 
nism should be observable in Er MOTs when larger 
MOT beam intensities are employed. The highly mag- 
netic atom Ho (9 fj,B, Ag/gg = 1.7%) — of interest for 
quantum information processing |23j — would likely ex- 
hibit anisotropic sub-Doppler cooling as well, if it were 
proven amenable to repumper-less magneto-optical trap- 
ping. 
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